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Abstract 

Background: Many insects are chennically defended against predatory vertebrates and invertebrates. Nevertlieless, 
our understanding of tlie evolution and diversity of insect defenses remains limited, since most studies have 
focused on visual signaling of defenses against birds, thereby implicitly underestimating the impact of insectivorous 
insects. In the larvae of sawflies in the family Tenthredinidae (Hymenoptera), which feed on various plants and 
show diverse lifestyles, two distinct defensive strategies are found: easy bleeding of deterrent hemolymph, and 
emission of volatiles by ventral glands. Here, we used phylogenetic information to identify phylogenetic 
correlations among various ecological and defensive traits in order to estimate the relative importance of avian 
versus invertebrate predation. 

Results: The mapping of 12 ecological and defensive traits on phylogenetic trees inferred from DNA sequences 
reveals the discrete distribution of easy bleeding that occurs, among others, in the genus Atholia and the tribe 
Phymatocerini. By contrast, occurrence of ventral glands is restricted to the monophyletic subfamily Nematinae, 
which are never easy bleeders. Both strategies are especially effective towards insectivorous insects such as ants, 
while only Nematinae species are frequently brightly colored and truly gregarious. Among ten tests of phylogenetic 
correlation between traits, only a few are significant. None of these involves morphological traits enhancing visual 
signals, but easy bleeding is associated with the absence of defensive body movements and with toxins occurring 
in the host plant. Easy bleeding functions through a combination of attributes, which is corroborated by an 
independent contrasts test indicating a statistically significant negative correlation between species-level 
integument mechanical resistance and hemolymph feeding deterrence against ants. 

Conclusions: Our analyses evidence a repeated occurrence of easy bleeding, and no phylogenetic correlation 
including specific visual signals is significant. We conclude that the evolution of chemically-based defenses in 
tenthredinids may have been driven by invertebrate as much as by avian predation. The clear-cut visual signaling 
often encountered in the Nematinae would be linked to differential trends of habitat use by prey and predators. 
Further studies on (prey) insect groups should include visual signals and other traits, as well as several groups 
of natural enemies, to better interpret their relative significance and to refine our understanding of insect 
chemical defenses. 
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Predatory birds and ants. Predator-prey interactions. Diversity, Evolution 
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Background 

Insects live under the Sword of Damocles, since numer- 
ous vertebrate and invertebrate predators attempt to eat 
them [1,2]. Predation is therefore a primary driving force 
in the evolution of insects, which survive biotic attacks 
among others by chemically based defense strategies, and 
an intriguing interspecific diversity in defense strategies 
is observed {e.g,, [3-5]). A specific defense strategy varies 
during ontogeny, and relates to an adapted phenology, 
behavior, morphology, physiology, and/or chemistry [6,7]. 
Defense strategies of living organisms are shaped by 
evolutionary conservatism and ecological factors, but 
few studies have attempted to estimate the relative im- 
portance of each of these two influences by a large-scale 
analysis of a given insect group [8-11]. This is under- 
standable, since 'eco-evo' processes of systems including 
insect prey and their predators are intrinsically complex 
[12]. We emphasize here three major points contribut- 
ing to this complexity. 

First, numerous insects are herbivorous, which gives 
them the possibility to reallocate toxic or harmful plant 
compounds to their own benefit (Figure 1). Sequestration is 
the uptake and accumulation of exogenous allelochemicals 
in specific organs [13], but other possible fates of plant 
allelochemicals are, for example, their detoxification or 
excretion by the insect [14]. Further, defense chemicals 
can be produced endogenously [15]; such de novo pro- 
duction can occur in non-herbivores, but surprisingly 
also in herbivores feeding on plants containing deleteri- 
ous allelochemicals. Species may benefit from this by 
becoming more independent from the plant, and by 
combining exo- and endogenous production, insects can 
facilitate their shifts to novel host-plant species [10,16,17]. 



Second, numerous insects prey on other insects, and 
such species exhibit fundamental differences in their hunt- 
ing strategy as compared to insectivorous vertebrates. 
Even though some predatory insects are visual hunters, 
most tend to locate and identify potential prey primarily 
by means of olfactory and gustatory cues [18,19]. This 
contrasts with vertebrate predators such as birds, which 
almost exclusively rely on vision when foraging [20-23], 
even if tasting is an important second step [24]. The 
point is that we perceive our environment as birds do, 
prevalently by sight, which may explain why many studies 
focus on visual signals such as crypsis, aposematism 
and its often associated traits, gregariousness and mim- 
icry. Thus, ecological factors determining the evolution 
of chemical defenses in insects are less studied than the 
signaling of such defenses [25] (Figure 1). 

Third, defensive chemicals are often multifunctional. 
Bioactive compounds can be general irritants acting on 
the peripheral sensory system, or toxins of specific 
physiological action [26]. Chemically, they roughly cor- 
respond to volatiles and water-soluble compounds, re- 
spectively. An advantage (for the emitter) of volatiles is 
that they keep the predator at a distance, whereas the 
action of water-soluble compounds requires ingestion 
or at least contact by the predator; repellence is defined 
here as involving the olfactory system, whereas feeding 
deterrence the gustatory one [27]. However, all such 
chemical and functional distinctions remain quite arbi- 
trary. Defensive chemicals in one species are often a 
mixture of chemicals and can be multifunctional by 
including chemical precursors, solvents, and/or wetting 
agents of the active compounds, by showing a feeding 
deterrence and toxicity, or a repellent and topical activity. 
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Figure 1 Evolutionary interactions among trophic levels influencing chemical defensive strategies in phytophagous insects. 

Phytophagous insects are held in 'ecological pincers' consisting of top-down as well as bottom-up selective pressures in the case of host plants 
containing deleterious chemicals {red arrows). However, the insects may sequester plant compounds, and/or produce defensive chemicals 
themselves, and they can also combine chemical with non-chemical defensive traits, which are all traits eventually used upon attack by natural 
enemies (green arrows). 
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etc. [4,5,15,28-31]. Even a single compound can be 
multifunctional [32], and different compounds often act 
in synergy [33]. More generally, dose-dependent effects 
of a chemical are ubiquitous, as already observed about 
500 years ago by Paracelsus {e.g., [34-36]). Finally, the 
interspecific activity of allelochemicals have led to a 
subset of names and definitions depending on the 
beneficial/detrimental action of the compounds for the 
emitter versus receiver, but again, a given compound 
can fulfill several of such ecological functions [37]. 

To better understand the evolution of chemical defen- 
sive strategies in phytophagous insects, we aimed to 
reconstruct the phylogeny of the Tenthredinidae saw- 
flies, which constitute the major group of herbivorous 
Hymenoptera, and which show a large diversity in life 
histories. Tenthredinids exhibit high intimacy with their 
host plant since females lay their eggs into the plant tissue 
[11]. Their larvae generally live freely on plant leaves and 
are preyed upon by numerous vertebrate and invertebrate 
predators [38]. Two distinct chemical defensive strategies 
are known among tenthredinid larvae. On the one hand, 
species in the subfamily Nematinae possess eversible ven- 
tral glands, which emit a volatile secretion that is probably 
aimed primarily against predatory insects and secondarily 
towards birds [39]. On the other hand, some tenthredinid 
species, especially those belonging to the blennocampine 
tribe Phymatocerini, are characterized by being able of 
easy bleeding! which is a phenomenon so far unknown 
from other insects and that is different from reflex bleed- 
ing [40] . In species able of easy bleeding, the larval integu- 
ment readily disrupts under exogenous mechanical stress 
at any point of the body [40-42], and the oozing hemolymph 
that contains sequestered plant secondary metabolites 
[14,43-45] is strongly feeding deterrent to biting preda- 
tors such as ants and wasps [40,43,46]. Comparative 
bioassays and modeling of the integument surface struc- 
ture indicate that easy bleeders are more effectively 
defended against such invertebrate predators than against 
birds [41,47]. Besides ventral glands and easy bleeding, 
alternative or complementary larval defenses include a 
developed pubescence, an integumental secretion layer 
[48,49], and an endophytic lifestyle by galling, rolling, 
mining or boring in different plant tissues [50,51]. More- 
over, there is diversity in the cryptic or aposematic ap- 
pearance, and level of gregariousness among tenthredinid 
larvae [39,52,53]. Such a large and diversified range of 
defensive devices within this insect group prompted us 
to search for evolutionary patterns, by seeking an ex- 
planatory framework of ecological factors that would 
account for this diversity. Therefore, we mapped eco- 
logical and defensive traits on phylogenetic trees, and 
tested correlations between character pairs, with the 
aim to infer the relative impact of invertebrates versus 
vertebrates in the evolution of chemically-based defenses. 



Our general hypothesis was that if vertebrates would be 
the main driver in this evolution, then at least some spe- 
cific visual signals should be correlated, at a large phylo- 
genetic scale, with an underlying chemical defense (see 
Figure 1). 

Methods 

Study group and taxon sampling 

Tenthredinidae is the largest sawfly family with over 
5,500 species described worldwide, covering all conti- 
nents except Antarctica [54]. Most tenthredinid species 
are dietary specialists: larvae typically feed on one or a 
few related plant genera [55]. The majority of the hosts 
belong to di- and monocotyledonous angiosperms, but 
numerous species feed on gymnosperms, ferns, horse- 
tails, and even mosses [55]. Generally, tenthredinids 
have been subdivided in seven, but more recently in the 
six subfamilies Allantinae, Blennocampinae, Heterarthrinae, 
Nematinae, Selandriinae, and Tenthredininae [54]. 

The sawfly species used in this study comprise 106 
tenthredinid species (Additional file 1), with representa- 
tives from all subfamilies and 29 of their major tribes. 
We primarily focused on species for which data on chem- 
ical defense traits but no accurate phylogenetic analyses 
were available, which is especially the case for the 
Phymatocerini, represented here by 22 exemplars. From 
10 non-tenthredinid sawfly families, 13 species were 
included in as outgroups. 

Sawfly adults were identified following Benson [56] and 
Zhelochovtsev & Zinovjev [57], and larvae with Lorenz 
& Kraus [48]. Specimens were stored in 100% ethanol 
at -20 or -80°C, and vouchers are kept at the Royal 
Belgian Institute of Natural Sciences (JLB coUection; 
Additional file 1). 

DNA extraction, PGR amplification and sequencing 

Total genomic DNA was extracted from legs or abdo- 
mens of adult sawflies or parts of larvae following a 
standard CTAB protocol. We amplified and sequenced 
two mitochondrial genes. Cytochrome b (Cytb) and 
Cytochrome oxidase I (COI), and the nuclear 28S ribo- 
somal gene. Cytb was amplified (and the PCR products 
sequenced in both directions) using primers CB- J- 10933 
and CB-N-11367 [58]. For Col, we used a modified ver- 
sion of the Simon et al. [58] primer Cl-J-1718 (5'-GGA 
GGA TTT GGA AAT TGA TTA ATT CC-3') in com- 
bination with the reverse primer mod-A2590 (5'-ACT 
GCT CCT ATT GAT AAT ACA TAA TG-3'; GM, own 
design). For 28S, the primers 28SF2 (5'-CAC GAG CCG 
ATA GCG AAC AAG T-3'; GM, own design) and 
28SB2 (5'-CCA AGG CCT CTA ATC ATT CGC T-3'; 
GM, own design) were used. PCR reactions contained 
10 mM Tris-HCl, 50 mM KCl, 1.5 mM MgCl2, 50 (iM 
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of each dNTP, 0.4 \iM of each primer, and 0.026 units/ (il 
of Taq polymerase (Amersham Bioscience). 

The PGR programs consisted of an initial denaturation 
step at 94°C for 4 min, followed by 30 cycles of 94°C 
for 60s, annealing at 52°C (Cytb and COI) or 54°C 
(28S) for 60s, and extension at 72°C for 2 min. The cy- 
cles were followed by a final extension step at 72°C for 
10 min. PGR products were purified using either alcohol 
precipitation or the Amersham Bioscience GFX PGR 
and Gel Band Purification Kit. 

Gytb was sequenced on an Amersham ALF express 
automatic sequencer using Gy5 labelled primers and the 
Amersham sequencing kit, and with an annealing temper- 
ature of 50°G. GOI and 28S products were sequenced using 
ABI PRISM BigDye Terminator cycle sequencing kits and a 
Perkin Elmer ABI sequencer at the Vlaams Instituut voor 
Biotechnologie in Antwerp, Belgium. 

Sequences were assembled and checked using the 
base-calling software of the respective sequencers, and 
then aligned using GlustalX v. 1.81 [59]. The alignments 
were corrected by eye. Alignment was straightforward 
for the mitochondrial genes, which contained very few 
insertions or deletions (indels), and in which codons 
could serve as reference. By contrast, numerous indels 
were present in the 28 S sequences, although these 
tended to occur mainly between outgroup and ingroup 
taxa. The final dataset contains 397, 862, and 999 
aligned base pairs for Gytb, GOI, and 28S, respectively 
(2,258 bp in all) (Additional file 2). All sequences have 
been submitted to GenBank under accession numbers 
KF528387-KF528662, and the full dataset (as well as 
resultant trees) are also available in TreeBase at http:// 
purl.org/phylo/treebase/phylows/study/TB2:S14547. 

Phylogeny reconstruction 

To reduce the effects of missing data, the full sequence 
alignment was split into two separate datasets: "Dataset 
1" included all 13 outgroup taxa and the 40 tenthredinid 
species that had sequences of all three genes (see 
Figure 2). "Dataset 2" included only outgroups from 
non-blasticotomid Tenthredinoidea (4 spp. representing 
Argidae, Pergidae, and Diprionidae), and all 106 ingroup 
taxa (see Figure 3). 

Both datasets were analyzed using Bayesian phylo- 
genetic inference as implemented in MrBayes v. 3.1.2 
[60] and BEAST v. 1.5.2 [61]. Prior to the runs, best- 
fitting substitution models for GOI (TVM+I+G), Gytb 
(TVM+I+G), and 28S (GTR+I+G) were identified under 
the Akaike information criterion in jModelTest 2.1.3. [62]. 
Because the TVM model is not implemented in MrBayes 
V. 3.1.2, we used a separate, unlinked GTR+I+G model of 
substitution for each gene in all phylogenetic analyses. 

The MrBayes analysis of Dataset 1 implemented de- 
fault priors and included two independent runs of four 



incrementally heated chains (t = 0.2) that were run for 
6 million generations, while sampling trees from the 
current cold chain once every 100 generations. The 
first 10,001 trees sampled prior to chain stationarity were 
discarded as a burnin from each run, and a Bayesian con- 
sensus tree showing all compatible groupings was calcu- 
lated on the basis of the 100,000 trees that remained in 
the combined tree sample. Runs with Dataset 2 were 
otherwise similar, but each run included six chains with 
the temperature parameter set to 0.1, and the analysis 
was run for 10 million generations. After deleting a 
burnin of 30,001 trees from both runs, a consensus tree 
was calculated on the basis of the remaining 140,000 trees. 

Topologically unconstrained BEAST runs of Dataset 1 
employed an uncorrected relaxed lognormal clock model 
of rate variation among branches, a Yule prior on speci- 
ation, and default priors for other parameters. Two 
independent runs with automatic tuning of operators 
were run for 60 million generations, while sampling 
trees and parameters every 1,000 generations. After 
discarding 10,001 trees from both runs as a burnin, 
the tree files were combined using LogGombiner (part of 
the BEAST package). A maximum clade credibility (MGG) 
tree showing mean node heights was then calculated 
on the basis of the 100,000 post-stationarity trees in 
TreeAnnotator (part of the BEAST package). Settings 
in the runs involving Dataset 2 were mostly identical 
(Additional file 3), but we ran four independent ana- 
lyses from which trees were combined after a burnin of 
10,001 trees. The combined tree file was then thinned 
by resampling trees every 4,000 generations, and an MGG 
tree was calculated based on the remaining 50,000 post- 
stationarity trees. 

Character coding and reconstruction of ancestral states 

In order to infer the evolutionary history of traits related 
to niche use and defense within Tenthredinidae, we col- 
lected data on diet breadths, host-plant associations, 
and host features, as well as on larval ecology, behavior, 
morphology, and chemically-based defensive strategies 
in the species included in the phylogenetic trees. We 
then coded the data as unordered binary or multistate 
characters, and reconstructed ancestral states by single- 
rate (Mkl) maximum likelihood (ML) optimization in 
Mesquite v. 2.75 [63]. To accommodate phylogenetic 
uncertainty, reconstructions were performed across 1,000 
post-burnin trees obtained by resampling trees from the 
Dataset 2 BEAST tree sample at regular intervals in 
LogGombiner, and results were summarized across the 
nodes of the MGG tree (Figure 3) from the same analysis 
(Additional file 4). ML optimization in Mesquite is not 
possible if taxa with unknown or polymorphic states 
are present in the focal character, so such species were 
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Posterior probability (%) (BEAST) 
Posterior probability (%) (MrBayes) 




Blennocampinae / 
Heterarthrinae / 
Allantinae 



Tenth redininae 



TBPh Monophadnus sp.2 [42] 
TBPh Monophadnus sp.4 [86] 
TBPh Monophadnus sp.5 [140] 
TBPh Monophadnus monticola [85] 
TBPh Monophadnus sp.B [84] 
THFe Metallus lanceolatus [22] 
TBPh Eutomostethus luteiventris [80] 
TBPh Eutomostethus ephippium [79] 
TBWa Waldheimia carbonaria [250] 
TBWa Halidamia affinis [90] 
TBBI Monardis plana [1 34] 
TAEm Harpiphorus lepidus [246] 
TBTo Tomostethus nigritus [4] 
THFe Notofenusa sp. [234] 
THCa Endelomyia aethiops [92] 
TAEm Ametastegia carpini [152] 
TTMa Macrophya ribis [109] 
TTMa Macrophya alboannulata [54] 
TTMa Pachyprotasis antennata [1 1 6] 
TTMa Pachyprotasis rapae [55] 
TTTe Tenth redo scrophulariae [1 4] 
TTTp Tenthredopsis tarsata [1 22] 
TTTp Aglaostigma sp. [119] 
TTSi S/ob/a sfurAn// [52] 
TAAI Allantus rufocinctus [225] 
TAAI Allantus viennensis [1 58] 
TSSe Brachythops flavens [1 36] 
TSSe Selandria serva [105] 
TSAn Aneugmenus pad! [1 1 ] 
TSDo Dolerus vestigialis [144] 
TSDo Dolerus aeneus [57] 
TSSt Stromboceros delicatulus [12] 
TNHo Hoplocampa testudinea [192, 30] 
TNPr Pristiphora geniculata [33] 
TNSu Susana cupressi [26] 
TNCI Cladius compressicornis [34] 
TAAt /\f/7a//a rosae [39] 
TAAt /\f/7a//a /*erta [66] 
TAAt Athalia cordata [65] 
TAAt Athalia Scutellariae [40] 
DIPRIONIDAE Neodiprion sert/fer [129] 
DIPRIONIDAE Gilpinia hercyniae [128] 
ARGIDAE /\rge ustulata [25] 
PERGIDAE Giladeus tuxius [232] 
BLASTICOTOMIDAE Blasticotoma ff//cef/ [137] 
PAMPHILIIDAE Acantholyda erythrocephala [37] 
PAMPHILIIDAE Neurotoma saltuum [59] 
CEPHIDAE Calameuta filiformis [127] 
XYELIDAE Xye/a cun/a [62] 
CEPHIDAE Janus compressus [61] 
MEGALODONTESIDAE Megalodontes thor [^30] 
XIPHYDRIIDAE Xiphydria prolongata [124] 
ORUSSIDAE Orussus abietinus [183] 

Figure 2 Relaxed molecular-clock phylogeny of the Tenthredinidae and selected outgroup taxa. The BEAST MCC tree is based on analysis 
of Dataset 1, wliicli includes tenthredinids that have sequences from all three genes, and all outgroups. Numbers above branches are posterior 
probabilities (%) from the BEAST analysis, numbers below branches show corresponding values from the MrBayes run (clades not present in the 
MrBayes tree are indicated by hyphens). Grey shaded bars show 95% highest posterior density intervals for relative node ages for nodes with 
posterior probabilities over 50%. 



Athaliinae 



deleted from the matrix and associated trees prior to 
each reconstruction. 

Host plants were compiled from literature sources 
[55,64] and each sawfly was coded for its diet breadth 
(Figure 3). In the case of species for which reliable host- 
plant records were missing, diet breadth was coded as 
unknown, unless revealed by further laboratory testing 
with larvae from the same population. 



To evaluate the toxicity of plants on which sawflies 
feed, each plant genus (and species, as far as possible) 
was associated with the occurrence of toxins, by referring 
to standard works on the chemistry of plants [65-70] and 
to smaller and/or more recent works {e.g., [45]). A plant 
taxon was considered toxic if the leaves contain secondary 
metabolites from one or more of the following chemical 
classes: alkaloids (including steroid alkaloids and steroid 
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(See figure on previous page.) 

Figure 3 Relaxed molecular-clock phylogeny of the Tenthredinidae, and the distribution of various larval ecological and defensive 
traits within the group. Tine BEAST MCC tree is based on analysis of Dataset 2, wliicli includes all sequenced tenthredinids as well as 
representatives from three non-blasticotomid families in Tenthredinoidea. Posterior probabilities (%) resulting from analyses in BEAST and MrBayes 
are given above and below branches, respectively (clades not present in the MrBayes tree are indicated by hyphens). Grey shaded bars show the 
95% highest posterior density intervals for relative node ages for nodes with posterior probabilities exceeding 50%. Branch colors denote host 
plant classes of the sawfly species (see legend) and ancestral reconstructions based on maximum-likelihood optimization across 1,000 post-burnin 
trees (see Additional file 4A). In the table to the right of the tree, diet breadth, plant toxicity, and defensive traits (from left to right) are coded as 
shown in Table 1. (?) Unknown; (x) not applicable. 



saponins, which are closely related in terms of biosyn- 
thesis and metabolism; [71]), glucosinolates, cyanogenic 
glycosides, and non-protein amino acids. All these com- 
pounds, except steroid saponins, have in common the 
inclusion of one or more nitrogen atoms. Nitrogen- 
containing secondary metabolites show acute toxicity 
and/or strong feeding deterrence towards vertebrates 
and/or invertebrates, and they are the most common 
defensive chemicals of plants [65,68-70,72]. A plant 
taxon was considered non-toxic if it only contains 
secondary metabolites that do not contain nitrogen, such 
as phenolics {e.g., coumarins, phenolic glycosides, and the 
widespread flavonoids), terpenoids {e.g., iridoid glycosides, 
triterpenoid saponins), or ranunculin (characteristic of 
the Ranunculaceae). Following the specific host plant(s) 
of each sawfly species, host toxicity was then coded as 
never' (code '0'), sometimes ('1'), or always' ('2'), depend- 
ing on the possible occurrence of toxins in the diet. For 
instance, the code was '0' for a specialist sawfly species 
feeding on a non-toxic plant genus, 1' for a generalist 
feeding on both toxic and non-toxic hosts, and '2' for a 
sawfly species only feeding on a toxic plant, or feeding on 
several plant taxa which are all toxic. 



Ten ecological traits linked to the behavior, morph- 
ology and chemical ecology of the sawfly larvae were 
coded as far as these traits are involved in defense (see 
Figure 3). The data were extracted from standard works 
on sawflies {e.g., [48,55,64,73] and literature therein), a 
specific work on easy bleeding [40], as well as unpub- 
lished observations and sources. For traits changing dur- 
ing successive larval stages, the last stage preceding the 
(often non-feeding) eonymph was considered. 

Correlation analyses 

The existence of phylogenetic correlations among various 
ecological and defensive traits was evaluated by Bayesian 
stochastic character mapping [74,75] as implemented in 
SIMMAP V. 1.5.2 [76]. For these analyses, we selected 10 
out of the 66 character-pair comparisons that are possible 
among the 12 focal traits listed in Table 1. Most correla- 
tions to be performed were selected based on previously 
proposed hypotheses (see [39,40,47] and Table 2). State- 
by-state associations between characters were evaluated 
based on the dij statistic, which measures co-occurrence 
of states / and ; across branches in relation to the 
expectation under independent evolution [75]. Overall 



Table 1 Plant features plus ecological and defensive traits of tenthredinid sawfly larvae used in reconstructing 
ancestral states and analyzing phylogenetic correlations 

Character (Code) state 



Diet breadth 


(0) one plant species or genus, (1) at least two plant genera but of one family, 




(2) plant genera of at least two families 


Plant toxicity 


(0) never, (1) sometimes, (2) always 


Mechanical plant protection 


(0) free-living larva, (1) leaf miner, (2) borer, (3) galler 


Placement on leaf 


(0) leaf edge, (1) leaf upper- and/or underside 


Gregariousness 


(0) solitary, (1) aggregated, i.e., larvae distributed on a plant, generally < 3 per leaf. 




(2) truly gregarious, i.e., larvae on one leaf or several adjacent leaves 


Defensive body movements 


(0) dropping easily and/or violent movements, (1) no, (2) raising abdomen 


Predominant body coloration 


(0) green, (1) white ventrally and green dorsally, (2) white or yellow, (3) brown-grey to black. 




or white ventrally and dark dorsally 


Distinct dark to black spots 


(0) absent, (1) present 


Exocrine ventral glands 


(0) absent, (1) present 


Body setation and protrusions 


(0) with very short setae and without long protrusions, (1) with setae > 1/6 as long as body diameter. 




(2) with protrusions or spines > 1/6 as long as body diameter 


Integumental wax layer 


(0) no, (1) yes 


Easy bleeding 


(0) no, (1) yes 
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Table 2 Overall phylogenetic correlations between various ecological and defensive characters (D) and associated P-values, 



estimated by Bayesian stochastic mapping across a sample of 500 post-burnin trees 


Ref. 


Character (code) 


Character (code) 


D 


P 


[40] 


Diet breadth (1) 


Plant toxicity (2) 


0.196 


0.010 




Plant toxicity (2) 


Mechanical plant protection (3) 


0.104 


0.056 


[40] 


Plant toxicity (2) 


Easy bleeding (12) 


0.260 


0.000 




Placement on leaf (4) 


Integumental wax layer (1 1) 


0.032 


0.198 


[39] 


Gregariousness (5) 


Defensive body movements (6) 


0.230 


0.000 


[39] 


Gregariousness (5) 


Dark spots on body (8) 


0.061 


0.468 


[39] 


Defensive body movements (6) 


Dark spots on body (8) 


0.078 


0.164 


[47] 


Defensive body nnovennents (6) 


Easy bleeding (12) 


0.444 


0.004 




Predominant body coloration (7) 


Body setation and protrusions (10) 


0.113 


0.048 


[40] 


Predominont body coloration (7) 


Easy bleeding (12) 


0.109 


0.024 



Most comparisons were selected by following the mentioned references (Ref.). Results from character-by-character analyses performed using a reduced dataset 
are given in italics, and P-values remaining statistically significant (at P<0.05) after Holm's sequential Bonferroni correction are given in bold. See Additional file 5 
for detailed results. 



character correlations were measured using statistic A 
which is the sum of the absolute values of individual di/s 
between characters [75] . 

Before the main analyses, parameter priors were deter- 
mined based on MCMC analyses following the approach 
of Schultz and Churchill [77] in SIMMAP. For binary 
characters, best-fitting priors for the bias parameter 
(= beta distribution parameter a) and the overall evo- 
lutionary rate parameters (= gamma distribution parame- 
ters a and p), were determined based on the Dataset 2 
BEAST MCC tree (Figure 3), which had been rescaled to 
a length of 1. MCMC runs were performed using the 
default number of distribution discretization categories 
(31 for the bias parameter and 60 for the rate parame- 
ters), cycles (100,000), sampling frequency (200), burnin 
(10,000), and upper rate bound (1,000). Results were 
extracted in R v. 2.14.0 [78] using the sumprmcmc.r 
script provided in the SIMMAP installation package. In 
the case of multistate characters, we used an empirical 
prior for the bias parameter, while rate-parameter priors 
were determined as described for binary characters. All 
characters were treated as unordered. 

In order to accommodate phylogenetic uncertainty, 
all correlation analyses were performed across 500 post- 
burnin trees (rescaled to a length of 1) obtained by 
regularly thinning the original Bayesian tree sample 
from the Dataset 2 BEAST runs in LogCombiner. Pair- 
wise correlation analyses were configured using the 
aforementioned bias and rate priors for each character. 
The number of samples, prior draws, and predictive 
samples were set to 1, meaning that both observed and 
predictive sample sizes equaled 500 for each character 
pair. For the three pairwise tests involving character 12 
(= easy bleeding), the high amount of missing data in 
this trait led to exceedingly long run times. Therefore, 
these analyses were instead based on a reduced dataset. 



in which the tree used for prior determination, as well 
as the ones in the 500-tree sample, were pruned to in- 
clude only those taxa that had a known state for easy 
bleeding. 

Independent contrasts test 

In a survey of nine species in the blennocampine tribe 
Phymatocerini, Boeve & Schaffner [40] found a statisti- 
cally highly significant negative correlation between in- 
tegument resistance and hemolymph deterrence, that is, 
the less the integument is resistant to a standardized 
mechanical stress, the more ants are deterred by the 
hemolymph, and vice versa. However, their analysis 
treated species as independent data points, which can 
potentially lead to spurious results due to phylogenetic 
non-independence of species [79]. Hence, we re-analyzed 
an expanded dataset using Felsensteins [80] independent 
contrasts method implemented in the PDAPiPDTREE 
package v. 1.15 [81] in Mesquite. These expanded analyses 
were based on 21 tenthredinid species for which both 
integument resistance and hemolymph deterrence had 
been measured [40]. The tree used in these analyses 
(a reduced version of the one shown in Figure 4A) was 
obtained by pruning the BEAST MCC tree in Figure 3. 

Results 

Phylogenetic trees 

The trees from the sequence data reveal evidence for 
the monophyly of the Tenthredinidae (Figures 2 and 3), 
as indicated earlier [82]. Outside Tenthredinoidea, how- 
ever, missing data in some outgroup representatives lead 
to clearly wrong groupings in Dataset 1 analyses, so the 
basal parts of the tree (Figure 2) should be treated with 
caution. This particularly concerns the placement of 
Xyelidae within Cephidae in the BEAST MCC tree, as 
well as the apparent polyphyly of the Pamphilioidea 
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Figure 4 Part of the phylogenetic tree of tenthredinids with estimated levels of traits linked to easy bleeding, and plot of independent 
contrasts extracted from a phylogeny that includes only species with no missing data. The tree in (A) was obtained by pruning tlie BEAST 
MCC tree in Figure 3, plots on tlie riglit-liand side of tine tree sliow levels of integument resistance and hemolymph deterrence estimated for the 
included species ([40,41] and U. Schaffner, unpublished data). Species excluded from the independent contrasts test due to missing data are 
denoted by gray terminal branches and parenthesized names. The scatterplot in (B) shows standardized contrasts for 21 nodes on the tree that 
include only species that have estimates for both traits, as well as the regression line forced through the origin. 



(Pamphiliidae + Megalodontesidae) in both analyses {cf, e.g., 
[83]). Within Tenthredinidae, the tree topologies are 
congruent in the monophyly and basal positioning of 
the genus Athalia, which justifies its placement in a dis- 
tinct subfamily, the Athaliinae, as proposed earlier {e.g., 
[55]). The topologies are also congruent in confirming 
the monophyly of the Nematinae [82]. Representatives 
of the Selandriinae, with the exception of the tribe 
Heptamelini, are grouped together. Each of the remaining 
traditional subfamilies, i.e., the AUantinae (with the 
aforementioned exclusion of Athalia), Blennocampinae, 
Heterarthrinae, and Tenthredininae, come out as poly- 
phyletic, and the groups are generally supported by low 
posterior probabilities. In traditional classifications, the 
AUantinae was, indeed, recognized very soon as an arbi- 
trary group [84], which is less the case for the three 
other subfamilies. However, in all subfamilies higher 
probabilities are obtained at lower-level (younger) clades, 
which allows the following conclusions. 

Easy bleeding is particularly common among a 
Blennocampinae tribe, the Phymatocerini ([40], Figure 3), 
which is a group defined by a combination of morpho- 
logical characters [73]. Our analysis does not demon- 
strate its monophyly (Figure 3) and rather shows two 
distantly related clades, one centered' on Monophadnus, 
and another on Rhadinoceraea. The latter clade includes 
Phymatocera and Paracharactus, and Eutomostethus is 
close to it. The wealdy supported exclusion of Monophadnus 
spinolae from other Monophadnus species, as well as the 



strong support for the grouping of Rhadinoceraea + 
Phymatocera + Par achar actus, are both reflected by mor- 
phological characters ([73], SMB, personal observation). 
The fact that the Phymatocerini are unique among the 
Blennocampinae in commonly feeding on plants containing 
steroidal saponins and alkaloids [40], which is clearly not a 
trait considered in the traditional classification of sawflies, 
lends additional support to the hypothetical monophyly of 
this tribe. 

Defense diversity 

A large diversity of lifestyles and defensive traits is 
found in tenthredinid larvae (Figure 3). Some traits 
evolved repeatedly, in at least two species groups, such 
as easy bleeding in Athalia and the Phymatocerini, leaf 
mining in the (possibly polyphyletic) Heterarthrini and 
Pseudodineurini, and an integumental wax layer in 
some Blennocampinae and Tenthredininae, and AUantinae 
(Additional fUe 4). In contrast, other traits are known 
from only one taxon. Examples are the eversible ventral 
glands in the Nematinae, the slimy covering in Caliroa, 
hemolymph spitting in Siobla, and fruit boring in 
Hoplocampa (Additional file 4). Moreover, a single species 
can combine at least two traits, for instance, aposematism 
and gregariousness, crypsis and a soUtary lifestyle, the 
presence of ventral glands and an endophytic lifestyle, 
or ventral glands and aposematism. However, easy bleed- 
ing and the presence of ventral glands never co-occur, 
meaning that no easy bleeder possesses ventral glands. 
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and that nematine species are never easy bleeders 
(Figure 3). The trees also indicate that easy bleeding 
appeared (and was lost) at least five times: in the 
Athaliinae, Allantinae, Selandriinae, Tenthredininae, and 
Blennocampinae (Phymatocerini), with a radiation of 
the phenomenon within the last of these taxa (Figure 3, 
Additional file 4). 

The wide range in overall diet breadth of tenthredinids 
impedes the recognition of a clear host- affiliation pattern 
for sawfly subgroups on host plant families and even 
orders. Most tenthredinid species feed on eudicots, 
with the two major exceptions that most Selandriinae 
feed on pteridophytes or monocots, and part of the 
Phymatocerini feed on monocots (Additional file 4). Plants 
containing toxic secondary metabolites are the host for 
species of Athalia, Selandriinae, (leaf-mining) Nematinae 
as well as the two Phymatocerini, Monophadnus- and 
Rhadinocemea'Centeiced, clades (Figure 3, Additional file 4). 

Associations among traits 

From the ten chosen pairwise comparisons, six yielded 
statistically significant overall correlations, but only three 
of them remain significant after Holms sequential 
Bonferroni correction: plant toxicity with easy bleeding, 
gregariousness with defensive body movements, and such 
movements with easy bleeding (Table 2, Additional file 5). 
More specifically, the results indicate that plant toxicity 
is associated with easy bleeding, easy bleeding with the 
absence of defensive body movements, a solitary habit 
with dropping and/or violent movements, aggregation 
with the absence of defensive movements, and true gre- 
gariousness with raising abdomen (Additional file 5). 

Felsenstein s independent contrasts test revealed a statis- 
tically significant negative correlation between species- 
level integument resistance and the rate of hemolymph 
deterrence (r = -0.393, = 0.155, P = 0.039; Figure 4B). 

Discussion 

The description and analysis of chemical defense mecha- 
nisms across insects, mainly in lepidopteran and coleop- 
teran herbivores, initiated the search for general trends 
in the taxonomic distribution and evolution of such 
mechanisms. Research using empirical and manipulative 
tests on predator-prey systems, computational modeling, 
and phylogeny-based approaches has identified sequential 
steps in the evolution of prey defensive traits as well as 
plant-insect interactions {e.g., [8,14,85-90]). However, nearly 
all such studies, even when they embrace multitrophic 
interactions at once, focus explicitly or implicitly on 
(dis) advantages as well as evolutionary sequences and 
consequences of visual prey signals. In this context, 
there is good evidence that the evolution of aposemat- 
ism is accompanied by an increased diversification of 
lineages, as shown by paired sister-group comparisons 



in insects and other animal taxa [91]. Further, chemical 
adaptation (unpalatability) preceded morphological (warn- 
ing coloration) and behavioral (gregariousness) adapta- 
tions in insects [8,85,87,89,92]. However, the next step in 
understanding the evolution and diversity of insect chem- 
ical defenses is to explain how unpalatability itself evolved, 
which remains a largely unexplored question. 

Since distastefulness in aposematic phytophagous in- 
sects often relies on plant chemistry, dietary specialization 
would favor aposematism due to physiological processes 
needed to cope with the ingested toxins [14,93]. Chemical 
specialization that is not necessarily related to plants' 
taxonomic affiliation also promotes aposematism, while 
similar chemical profiles of secondary compounds across 
plant taxa facilitate niche shifts by phytophagous insects 
[10,93,94], which in turn may enhance the diversity of 
chemicals underlying aposematism. But, shifts in resource 
or habitat are probably less common than previously 
expected, as shown for sawfly larvae and caterpillars [95,96], 
and all aforementioned considerations are true for exogen- 
ous but not endogenous insect toxins, because these are 
per se unrelated to host affiliation. By the examination 
of an insect group with defensive features including, 
among others, bright and cryptic colorations, we could 
evaluate the probably literature-biased prevalence of 
avian over invertebrate predation in the evolution of 
insect defensive strategies. 

Our study reveals a diversity of defensive strategies. 
The tree-based analyses confirm previous conclusions 
from chemical and bioassay analyses on selected plant- 
insect systems. First, easy bleeding is a defense based on 
toxins from plants [40,43-45] since easy bleeders tend to 
feed on plants that contain such chemicals (Additional 
file 5). Second, easy bleeders move slowly and become 
immobfle once disturbed, whereas other defensive be- 
haviors are associated with non-easy bleeders ([47], 
Additional file 5). A third significant association to arise 
from our analyses is between defensive body movements 
and gregariousness. These behaviors are components of 
visual signals, but they impact predator-prey interac- 
tions also physically. They are effective towards birds as 
well as invertebrates, for instance, when attacking ants 
can be knocked or dislodged by defensive body move- 
ments, or when foraging ants ignore the presence of an 
immobile larva [39,47,97]. Associations including more 
specifically visual traits of the tenthredinid larvae were 
expected to be significant. Each of the two traits, 'dark 
spots on body and predominant body coloration', was 
tested against two other traits, but, surprisingly, none 
of these four associations is significant (Table 2). In 
particular, 'dark spots on body', which contributes to a 
conspicuous coloration is not associated with gregari- 
ousness. This contradicts with studies on several insect 
groups, including the Nematinae, that emphasize the 



Boeve et al. BMC Evolutionary Biology 201 3, 13:1 98 
http://www.bionnedcentral.conn/1471 -21 48/1 3/1 98 



Page 11 of 14 



link between aposematic coloration and gregariousness 
[8,9,39,52,87,98]. It seems that by studying the vast group 
of the tenthredinids we incorporated multiple defensive 
traits related to visual, mechano-physical and/or chemical 
cues, without focusing only on those known a priori to 
be directed against birds. 

Besides insectivorous birds, predatory insects and 
especially ants are known to shape communities and 
influence the evolution of sawflies [99,100], and besides 
easy bleeding, a second main defensive strategy is the 
presence of volatile-emitting ventral glands. Both easy 
bleeding and ventral glands are most effective against 
predatory insects such as ants, and less so against birds 
[39,47]. We assume intertwined roles played by inverte- 
brate and vertebrate predators on the evolution of de- 
fensive strategies in tenthredinids. Their basal taxon, 
Athalia, and other taxa use easy bleeding as defense, 
and the tenthredinids has radiated into species-rich 
groups such as the Selandriinae (970 species), Nematinae 
(1,250), and Tenthredininae (1,720) [54], which illustrates 
the success of the family. Predation is generally believed 
to be a main driver in the evolution of insects, and the 
observed patterns suggest that the evolution and radiation 
of several tenthredinid subgroups have been driven by 
invertebrate rather than by vertebrate predators, and by 
which easy bleeding arose as a first defensive strategy. It 
remains unknown why this unique defensive strategy 
did not evolve in other insects while it was gained 
and lost several times in tenthredinids (Figure 3 and 
Additional file 4). Conversely, the use of a volatile secre- 
tion produced by exocrine glands is rather common in 
insects [4], but within the tenthredinids it is restricted 
to the Nematinae, and defense by ventral glands there- 
fore seems to be an alternative to easy bleeding. To be 
effective, the two defense strategies require quite oppos- 
ite behaviors, by raising the abdomen and by becoming 
immobile, respectively, which may explain why they 
evolved in a mutually exclusive way. In contrast, more 
exclusive visual signals could theoretically complement 
both of them. 

It is then intriguing that the Nematinae include rela- 
tively frequent cases of brightly colored and truly gre- 
garious species ([48,53], Figure 3), which indicates a 
more specific evolutionary impact of birds. Since both 
easy bleeding and ventral glands are primarily directed 
against invertebrate predators [39-41,47], the paradox is 
that only the latter defense is repeatedly linked to apo- 
sematic coloration, while the chemistry underlying both 
defenses is potentially as effective against invertebrates 
as against vertebrates (see later), and volatiles are not 
particularly effective against birds. Athalia, Nematinae, 
and also Phymatocerini appeared within a relatively 
short time range or even concurrently (Figures 2 and 3), 
so that it is not likely that one defense strategy arose 



much later than the other, in response to a new predation 
type. The Nematinae compared to all other tenthredinids 
are however singular in having their greatest diversity in 
the northern Holarctic, and they have a propensity to 
feed on willows (Salix) and other trees and shrubs, 
whereas most Athalia feed on Brassicaceae, and most 
Phymatocerini on Ranunculales and Liliales [55,73,95], 
the three latter groups of plants being herbaceous. In- 
sectivorous birds forage commonly in open (understory 
and canopy) forest habitats and probably less often at 
ground-level where they themselves suffer a higher 
predation risk [101-103], whereas ants occur more 
equally across all vegetation levels [104]. Differential 
trends of microhabitat-linked predation pressures, by 
ants and birds on Nematinae versus more prominently 
by ants on Athalia and Phymatocerini, may have driven 
the evolution of differing defensive traits [105]. The 
underlying chemicals, water-soluble compounds and 
volatiles, are dissimilar, too. Both types of chemicals can 
act on invertebrates as well as vertebrates [4,31,68,72]. 
Since any predator logically approaches before it attacks 
a potential prey, volatiles by acting at distance are more 
prone than water-soluble compounds to prevent an at- 
tack, but the latter compounds have a more profound 
physiological effect upon ingestion [26]. The effective- 
ness of different types of allochemicals is moreover 
affected by the predator-prey body size ratio and the 
consequently possible set of behavioral interactions 
between both protagonists [106,107]. It is within one 
type of allelochemicals that the diversity of selective 
pressures imposed by predators may promote preys' 
chemical diversity, while the type itself of chemicals 
would be determined by basic, morphological and 
physiological features. How predators promote chemical 
defense diversity requires further analyses by focusing 
not only on birds [90] but also insectivorous insects. 
Tenthredinids are a singular group of prey insects due 
to the unique occurrence of easy bleeding, but our case 
study on them evidences general patterns of chemically- 
based prey adaptations, and it adds to our overall under- 
standing of chemical defense diversity in insects. 

Conclusions 

Contrasting selective pressures imposed by various 
natural enemies on insect herbivores are likely to lead 
to the evolution of distinct defensive syndromes that 
potentially can be identified based on phylogenetic cor- 
relations among multiple independent traits. In the 
family Tenthredinidae, a staggering diversity of defen- 
sive strategies has evolved, and our macro-evolutionary 
analyses uncover several cases of evolutionary non- 
independence among anti-predator traits. In the par- 
ticular case of easy bleeding, an independent contrast 
test confirmed the existence of a negative phylogenetic 
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correlation between the mechanical resistance of the 
integument and the hemolymphs feeding deterrence to- 
wards ants. Since water-soluble compounds from the 
hemolymph of easy bleeders {Athalia and Phymatocerini 
species, among others) as well as volatiles from the 
ventral glands (in the Nematinae) are more prone to 
act as a defense against predatory insects than birds, it 
is likely that the obvious visual signaling often encoun- 
tered in the Nematinae is caused indirectly by differen- 
tial trends of habitat use by sawfly prey versus predator 
groups. Although several ecological and defensive traits 
were screened in tenthredinid larvae, none of those 
referring to specific visual signals were significantly 
correlated with the others. We conclude that, without 
neglecting the selective pressure by insectivorous birds, 
it seems necessary to emphasize the overall evolutionary 
impact of invertebrate predators on insect defensive 
strategies. 

Several theories on plant-insect relationships account 
for the diversity of plant defenses [88], while extending 
such theories to predator-prey relationships is much 
rarer [86] because especially those interactions involving 
herbivorous prey can become extremely complex. The 
mechanisms of arms races between predators and danger- 
ous prey imply coevolution rather than escalation [108]. 
From a predator s perspective, coevolution and escalation 
differ in what selective agents are responsible: the 
defense of a prey, or the attack by a predator. From a 
tenthredinid s perspective, the present study and others 
indicate that sawfly species face guilds of vertebrate and 
invertebrate predators, but also parasitoids and patho- 
gens [109-113]. Chemicals conform but also differ in 
their bioactivity on distant taxonomic groups such as 
invertebrates and vertebrates. A partially variable bio- 
activity can 'bridge' the use of defensive chemicals from 
one target group to another, and, hence, promote chemical 
diversification. A possible pattern of the macroevolution of 
insect chemical defenses would be that allelochemicals 
effective on invertebrates were co-opted for their bio- 
activity on birds. There is a need for further research on 
such adaptive cascades in insects. 
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